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Abstract 

We study nuclear shadowing for transverse and longitudinal photons. The coherence 
length, which controls the onset of nuclear shadowing at small Bjorken longer 
for longitudinal than for transverse photons. The light-cone Green function technique 
properly treats the finite coherence length in all multiple scattering terms. This is 
especially important in the region xsj > 0.01, where most of the data exist. NMC 
data on shadowing in deep inelastic scattering are well reproduced in this approach. 
We also incorporate nonperturbative effects, in order to extrapolate this approach to 
small photon virtualities Q^, where perturbative QCD cannot be applied. This way, 
we achieve a description of shadowing that is based only on quark and gluon degrees 
of freedom, even at low Q^. 
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1 Introduction 



The use of nuclei instead of protons in high energy scattering experiments, such as deep 
inelastic scattering (DIS), provides unique possibilities to study the space-time development 
of strongly interacting systems. In experiments with proton targets the products of the 
scattering process can only be observed in a detector which is separated from the reaction 
point by a macroscopic distance. In contrast to this, the nuclear medium can serve as 
a detector located directly at the place where the microscopic interaction happens. As a 
consequence, with nuclei one can study coherence effects in QCD which are not accessible 
in DIS off protons nor in proton-proton scattering. 



Figure 1: At low xbj and in the target rest frame, the virtual photon 7* converts into a 
gg-pair long before the target. 

At high energies, nuclear scattering is governed by coherence effects which are most easily 
understood in the target rest frame. In the rest frame, DIS looks like pair creation from a 
virtual photon, see Fig. ^ Long before the target, the virtual photon splits into a gg-pair. 
The lifetime Ic of the fluctuation, which is also called coherence length, can be estimated 
with help of the uncertainty relation to be of order ~ 1/mj^XBji where xbj is Bjorken-x and 
mjv ~ 1 GeV is the nucleon mass. The coherence length can become much greater than the 
nuclear radius at low XBj- Multiple scattering within the lifetime of the gg fluctuation leads 
to the pronounced coherence effects observed in experiment. 

The most prominent example for a coherent interaction of more than one nucleon is the 
phenomenon of nuclear shadowing, i.e. the suppression of the nuclear structure function F2 
with respect to the proton structure function Fj', F2{xBj,Q'^)l{AFl{xBj,Q'^)) < 1, at low 
XBj ^0.1. Shadowing in low xbj DIS and at high photon virtualities is experimentally well 
studied by NMC p. 

What is the mechanism behind this suppression? If the coherence length is very long, 
as indicated in Fig. Q the gg-dipole undergoes multiple scatterings inside the nucleus. The 
physics of shadowing in DIS is most easily understood in a representation, in which the pair 
has a definite transverse size p. As a result of color transparency small pairs interact 

with a small cross section, while large pairs interact with a large cross section. The term 
"shadowing" can be taken literally in the target rest frame. Large pairs are absorbed by 
the nucleons at the surface which cast a shadow on the inner nucleons. The small pairs are 
not shadowed. They have equal chances to interact with any of the nucleons in the nucleus. 
From these simple arguments, one can already understand the two necessary conditions for 
shadowing. First, the hadronic fluctuation of the virtual photon has to interact with a large 
cross section and second, the coherence length has to be long enough to allow for multiple 
scattering. 
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2 Shadowing and diffraction in DIS 



Like shadowing in hadron-nucleus collisions, shadowing in DIS is also intimately related 
to diffraction The close connection between shadowing and diffraction becomes most 
transparent in the formula derived by Karmanov and Kondratyuk [Sj. In the double scat- 
tering approximation, the shadowing correction can be related to the diffraction dissociation 
spectrum, integrated over the mass. 
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da{YN XN) 
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Here 



dzpA{b, zY'l^^ 



is the formfactor of the nucleus, which depends on the coherence length 



(3) 



V is the energy of the 7* in the target rest frame and Mx is the mass of the diffractively 
excited state. The coherence length can be estimated from the uncertainty relation and is 
the lifetime of the diffractively excited state. If Zc 0, the shadowing correction in Eq. (|T|) 
vanishes and one is left with the single scattering term Ao'^ p. 



2.1 The dipole approach 

Note that Eq. is valid only in double scattering approximation. For heavy nuclei, however, 
higher order scattering terms will become important. These can be calculated, if one knows 
the eigenstates of the interaction. Fortunately, the eigenstates of the T matrix (restricted to 
diffractive processes) were identified a long time ago in QCD [3 E] as partonic configurations 
with fixed transverse separations in impact parameter space. For DIS, the lowest eigenstate is 
the qq Fock component of the photon. The total 7*-proton cross section is easily calculated, 
if one knows the cross section (Jqq{p) for scattering a gg-dipole of transverse size p off a 
proton. 




Here, a is the longitudinal momentum fraction carried by the quark in Fig. ^ The light- 
cone wavefunctions \l/^^^(a,p) describe the splitting of a transverse (T) and longitudinal (L) 
photon into a gg-pair. For small p, the light-cone wavefunctions can be calculated in pertur- 
bation theory (see e.g. [7j for explicit expressions), but at large p non-perturbative effects 
become important. In Ref. [H], these effects have been modeled by introducing a harmonic 
oscillator potential between the quark and the antiquark, which leads to a modification of the 
light-cone wavefunctions. The strength of this potential has been determined from data for 
photoabsorption on protons. This justifies the application of the dipole formulation at low 
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Q"^ and makes this approach an ahernative to the vector dominance model (see e.g. Ref. [H]). 
At large Q^, of course, the modified light-cone wavefunctions reduce to the perturbative ones. 

The dipole cross section is governed by nonperturbative effects and cannot be calculated 
from first principles. We use the phenomenological parameterization that is fitted to HERA 
data on the proton structure function. Note that higher Fock-states of the photon, containing 
gluons, lead to an energy dependence of a^g, which we do not write out explicitly. 

The diffractive cross section can also be expressed in terms of (Tqq. Since the cross section 
for diffraction is proportional to the square of the T-matrix element, |(7*|T|X)p, the dipole 
cross section also enters squared. 



/ 
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(5) 



where the brackets (. . .) indicate averaging over the light-cone wavefunctions like in Eq. (j3)). 
We point out that in order to reproduce the correct behavior of the diffractive cross section 
at large Mx, one has to include at least the qqG Fock-state of the 7*. This correction is, 
however, of minor importance in the region where shadowing data are available. 



2.2 The Green function technique 

If one attempts to calculate shadowing from Eq. with help of Eq. (0), one faces the 
problem that the nuclear form factor, Eq. (0), depends on the mass Mx of the diffractively 
produced state, which is undefined in impact parameter representation. Only in the limit 
Ic ^ Ra, where Ra is the nuclear radius, it is possible to resum the entire multiple scattering 
series in an eikonal-formula 

a^'^ = j (^1 - exp 

The nuclear thickness function T{b) = dz pA^b, z) is the integral of nuclear density over 
longitudinal coordinate z and depends on the impact parameter h. The condition 1^ ^ Ra 
makes sure that the p does not vary during propagation through the nucleus (Lorentz time 
dilation) and one can apply the eikonal approximation. 

The condition Ic ^ Ra is however not fulfilled in experiment. For the case Ic ~ Ra, 
one has to take the variation of p during propagation of the qq fiuctuation through the 
nucleus into account, see Fig. El A widely used recipe is to replace Mx Q^, so that 
Ic l/{2mNXBj) and one can apply the double scattering approximation. This recipe was, 
however, disfavored by our investigation ^U]- Moreover, there is no simple recipe to include 
a finite Ic into higher order scattering terms. 

In Ref. PU] (see also Ref. jTTj) a Green function technique was developed that provides the 
correct quantum-mechanical treatment of a finite coherence length in all multiple scattering 
terms. Like in Eq. the total cross section is represented in the form 

(jT*^ = Aa^'P - Act, (7) 
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Figure 2: Propagation of a gg-pair through a nucleus. Shown is the case of a finite coherence 
length, where the transverse motion is described by the Green function W {p2, zi). 

where Aa is the shadowing correction, 

00 00 
Aa = 2^^5Z / / dzipA{h,zi) j dz2pA{h,Z2) 



f da j d^p2 '^q^^{p2,a) aqq{p2)A{p2, zi, Z2, a 



with 



A{p2, zu Z2, a) = j d^pi W{p2, Z2\pi, zi) e-'""^""^''-'^^ aM ^^^(pi, a). (8) 
Here, 

1 Q'^ail - a) +m'^f 
1^'^'= 2ua{l-a) ^ ' 

is the minimal longitudinal momentum transfer when the photon splits into the qq dipole 
[mf is the quark mass). 

The shadowing term in Eq. ((71) is illustrated in Fig. |21 At the point Zi the photon 
diffractively produces the qq pair (7*iV — > qqN) with a transverse separation pi. The pair 
propagates through the nucleus along arbitrarily curved trajectories, which are summed 
over, and arrives at the point Z2 with a separation p2. The initial and the final separations 
are controlled by the light-cone wavefunctions a). While passing the nucleus the qq 

pair interacts with bound nucleons via the cross section Cgglp) which depends on the local 
separation p. The Green function W{p2, Z2\pi,zi) describes the propagation of the pair from 
Zi to Z2, see Eq. (jH}, including all multiple rescatterings and a finite coherence length. Note 
the diffraction dissociation {DD) amplitude, 

/dd(7* qq) = i^S^(Pi, a)aqq{pi), (10) 

in Eq. (jH|). At the position Z2, the result of the propagation is again projected onto the 
diffraction dissociation amplitude. The Green function includes that part of the phase shift 
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Figure 3: The xbj dependence of shadowing in DIS for the structure function of tin relative 
to carbon. The data are from NMC |Tj. The full curves are calculated in the Green-function 
approach, including the nonperturbative interaction between the q and the q. The dashed 
curve does not include this interaction. The figure is from Ref. [T^ . 



between the initial and the final photons which is due to transverse motion of the quarks, 
while the longitudinal motion is included in Eq. (jH)) through the exponential. 

The Green function W{p2, Z2', pi,zi) in Eq. (jH)) satisfies the two dimensional Schrodinger 
equation, 

. dW{p2,Z2\pi,Zi) ^A{p2) ^ f \ (u I- \ 

1 = - 7^ — 7--, r W[p2, Z2\pi, Zi) - - a{p2) pA{b, Z2) W{p2, Z2\pu Zl) (11) 

OZ2 Zva^i — a) Z 

with the boundary condition W{p2,zi\pi,zi) = 6^'^\p2 — pi). The Laplacian A(p2) acts 
on the coordinate p2- The kinetic term A/[2z/a(l — a)] in this Schrodinger equation takes 
care of the varying effective mass of the qq pair and provides the proper phase shift. The 
role of time is played by the longitudinal coordinate Z2. The imaginary part of the optical 
potential describes the rescattering. This equation has recently been solved numerically by 
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Nemchik [121 . 

The Green function method contains the Karmanov-Kondratyuk formula Eq. (0) and the 
eikonal approximation Eq. © as hmiting cases. In order to obtain the eikonal approximation, 
one has to take the hmit v oo. In this case, the kinetic energy term in Eq. (|TT|l can be 
neglected and with g™" — ^ one arrives after a short calculation at Eq. 0. One can also 
recover the Karmanov-Kondratyuk formula, when one neglects the imaginary potential in 
Eq. Then W becomes the Green function of a free motion. 

Calculations in the Green function approach are compared to NMC data in Fig. El Note 
that the leading twist contribution to shadowing is due to large dipole sizes, where nonper- 
turbative effects, such as an interaction between the q and the g, might become important 
(see above). Therefore, the solid curve is calculated with the modified light-cone wave- 
functions of Ref. [SI. The dashed curve is calculated with the conventional, perturbative 
light-cone wavefunctions, but including a constituent quark mass. Both curves are in rea- 
sonable agreement with the data. We stress that the Green function technique also takes 
into account some higher twist corrections to shadowing. This is essential for a successful 
description of the NMC data. In fact, it has been demonstrated in Ref. Jll that the leading 
twist approximation only poorly reproduces NMC data for shadowing in calcium. 

Note that for the data shown in Fig. El the coherence length is of order of the nuclear 
radius or smaller. Indeed, shadowing vanishes around Xsj ~ 0.1, because the coherence 
length becomes smaller than the mean internucleon spacing. Therefore, the eikonal approx- 
imation, Eq. 0, cannot be applied for the kinematics of NMC and a correct treatment of 
the coherence length becomes crucial. We emphasize that the calculation in Fig. 01 does not 
contain any free parameters. Following the spirit of Glauber theory, all free parameters are 
adjusted to DIS off protons. 

2.3 Higher Fock states and the leading twist gluon shadowing 

The shadowing for quarks discussed in the previous subsection is dominated by the trans- 
verse photon polarization. Longitudinal photons, on the other hand, can serve to measure 
the gluon density because they effectively couple to color-octet-octet dipoles. This can be 
understood in the following way: the light-cone wave function for the transition 72 — > qq does 
not allow for large, aligned jet configurations as is the case for transversely polarized photons. 
Thus, all qq dipoles from longitudinal photons have size and the double-scattering term 
vanishes oc l/Q^. The leading-twist contribution for the shadowing of longitudinal photons 
arises from rescattering of the \qqG) Fock state of the photon. Here again, the distance 
between the q and the q is of order l/Q^, but the gluon can propagate relatively far from 
the gg-pair. In addition, after radiation of the gluon, the pair is in an octet state. There- 
fore, the entire |ggG)-system appears as a GG-dipole, and the shadowing correction to the 
longitudinal cross section can be identified with gluon shadowing. 

A critical issue for determining the magnitude of gluon shadowing is the distance the 
gluon can propagate from the gg-pair in impact parameter space, i.e. knowing how large 
the GG dipole can become. This value has been extracted from single diffraction data in 
hadronic collisions in Ref. [SI because these data allow the diffractive gluon radiation (the 
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triple-Pomeron contribution in Regge phenomenology) to be unambiguously singled out. The 
diffraction cross section (oc p^) is even more sensitive to the dipole size than the total cross 
section (oc p^) and is therefore a sensitive probe of the mean transverse separation. It was 
found in Ref. |H] that the mean dipole size must be of the order of Tq = 0.3 fm, considerably 
smaller than a light hadron. A rather small gluon cloud of this size surrounding the valence 
quarks is the only way that is known to resolve the long-standing problem of the small size 
of the triple-Pomeron coupling. The smallness of the GG dipole is incorporated into the LC 
approach by a nonperturbative interaction between the gluons. 

We calculate gluon shadowing as function of x at fixed and as a function of at 
fixed x, integrated over the impact parameter h. The results are shown in Fig. EJ In the 
left-hand plot, one observes that gluon shadowing vanishes for x > 0.01. This happens 
because the lifetime of the -fluctuation becomes smaller than the mean internucleon 
distance of ~ 2 fm as a; exceeds 0.01. Indeed, in Ref. ^3] an average coherence length of 
slightly less than 2 fm was found for the |ggG)-state ai x = 0.01 and large ^ ^/'^o- 
Note that gluon shadowing sets in at a smaller value of x than quark shadowing because 
the mass of a -state is larger than the mass of a -state. We also point out that 
gluon shadowing is weaker than quark shadowing in the x-range plotted, because the small 
size of the GG-dipole overcompensates the Casimir factor in the GG-proton cross section, 
c"GG = (9/4)(Tgg. Note however that at this time, almost nothing is known about gluon 
shadowing from experiment, and theoretical approaches differ vastly (see e.g. Ref. jTHj for 
a comparison of different models). The plot on the right-hand side of Fig. |^ shows the Q^- 
dependence of gluon shadowing and clearly demonstrates that gluon shadowing is a leading- 
twist effect, with Rq only very slowly (logarithmically) approaching unity as — > oo. 

2.4 The mean coherence length 

The importance of the coherence length 1^ was already mentioned above: this quantity con- 
trols the onset of quark and gluon shadowing at small x bj and therefore governs shadowing 
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in the kinematical region that is of interest for most experiments. The Green function tech- 
nique is the only known way to include the quantum mechanically correct coherence length 
into all multiple scattering terms. However, since Ic is not well defined in impact parameter 
space, it appears only implicitly in the Green function approach. It is therefore useful to 
introduce the concept of a mean coherence length, as it was done in Ref. jl5j . 

A photon of virtuality and energy u can develop a hadronic fluctuation for a lifetime, 

2u P 



Ir 



(52 + XBj rriN 
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where Mgg is the effective mass of the fluctuation, and the factor P ^ 
usual approximation is to assume that M^^ ^ 



{1 + My Q'). The 
since is the only large dimensional scale 

available. In this case P = 1/2. 

The effective mass of a non-interacting gg-pair is well defined, M^^ = [m'j + p'^) / a{l — a), 
where Pt and a are the transverse momentum and fraction of the light-cone momentum of 
the photon carried by the quark, respectively. Therefore, P has a simple form, 

g2a(i 



a 



(13) 



where = a{l — a)Q'^ + mj. 

To find the mean lifetime of those fluctuations that contribute to shadowing, one should 
define the averaging procedure as 





P{Pt, a) 


fir qq)) 


{fir - 


> qq) 


fir 


qq)) 



(14) 



where fir qq) is the amplitude of diffractive dissociation of the virtual photon on a 
nucleon, Eq. (fTIUl . That way, P is weighted with the interaction cross section squared (T^gip) 
in the averaging procedure. Then, the mean value of the factor P{a,pT) reads for transverse 
and longitudinal photons. 



Jq da J (Ppid?p2 




CTqq (P2) P ip2 - Pi, «) ^J"^^ (Pl, «) (^qq (Pl) 


lo da J (Pp 


(pa) CTqq (P) 
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^(P2 -Pi,a) 



(PpT_ 

{2n)' 



exp i-ipT ■ ip2 - Pi))P («, p). 



(15) 



(16) 



Numerical results are shown in Fig. El Note that the nonperturbative interaction also 
modifies the invariant mass of the qq pair from Eq. (fT!^ . This was taken into account in 
Ref. ^H]- One observes that the mean coherence length for longitudinal photons is ap- 
proximately twice as long as for transverse photons. However, shadowing for the qq Fock 
component of a longitudinally polarized photon is higher twist (see above). The coherence 
length for the qqG Fock state of a longitudinal photon, which gives rise to leading twist 
gluon shadowing, is much shorter, resulting in a delayed onset of shadowing for gluons. 
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Figure 5: The mean coherence length as function of Q^. The curves are for qq fluctuations 
of transverse and longitudinal photons, and for qqG fluctuation, from the top to bottom, 
respectively. Dotted curves correspond to calculations with perturbative wavefunctions and 
an approximate dipole cross section oc p^. Dashed curves are the same, except the realistic 
parameterization of o-qq{p) is employed (see Ref. [T3] for details). The solid curves are 
calculated with the nonperturbative wavefunctions of Ref. jH] 

3 Summary 

DIS at low XBj is most naturally described in the color dipole formulation, because partonic 
configurations with fixed transverse separations in impact parameter space are eigenstates of 
the interaction. This salient feature allows one to calculate multiple scattering effects, such 
as nuclear shadowing, in a very easy way: at very high energies, one can simply eikonalize the 
dipole cross section. At realistic energies, however, corrections due to the finite lifetime Ic of 
the gg-pair become important. In Ref. jTUj, we succeeded in generalizing the Gauber-Gribov 
theory of nuclear shadowing by incorporating a finite Ic into all multiple scattering terms, 
using a Green function technique. Since nuclear shadowing is dominated by large dipole 
sizes, a nonperturbative interaction in form of an harmonic oscillator potential between the 
quarks was introduced in Ref. jH]. This makes the dipole approach applicable at low Q^. 
Note however that it is not necessary for this model to reproduce the vector meson masses 
or the coherence length of the vector dominance model |15j . 

The main nonperturbative input to all formulae, the dipole cross section, cannot be cal- 
culated from first principles. Instead we use a phenomenological parameterization for this 
quantity, which is determined from low XBj DIS. In the spirit of Glauber theory, nuclear 
effects are then calculated without introducing any new parameters. This way a good de- 
scription of NMC data on shadowing in DIS is achieved. We did not attempt to include 
antishadowing, since this effect probably is beyond the standard shadowing dynamics. 

The parameterization of the dipole cross section effectively also includes higher Fock 
states (containing gluons) through its energy dependence. These higher Fock states are 
however excluded from nuclear effects, if one eikonalizes only the qq cross section. The 
most striking consequence is that gluon shadowing [i.e. shadowing for longitudinal photons) 
appears to be higher twist. This problem is overcome by calculating the rescattering of the 
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|ggG)-Fock state of the virtual photon, which gives rise to the leading twist gluon shadowing. 
A detailed calculation is published in Refs. jHlEl- We emphasize that gluon shadowing sets 
in at smaller xbj than shadowing for quarks, because the larger mass of the |ggG')-Fock state 
leads to a shorter coherence length. This result is supported by a calculation of the mean 
coherence length ^H]- Shadowing disappears, when the coherence length becomes shorter 
than the mean internucleon separation (~ 2fm). 

The main advantage of the dipole formulation and our motivation for pursuing this 
approach is the insight it provides into the dynamical origin of nuclear effects, which are 
calculated without free parameters. Note that a variety of other processes can be described 
in the dipole language, e.g. Drell-Yan (DY) dilepton production [THI, gluon radiation [Til] 
heavy quark (201 and quarkonium production |25- A parameter free calculation of nuclear 
shadowing for DY for example is needed, if one aims at extracting the energy loss of a fast 
quark propagating through nuclear from DY data [22] • Furthermore, this approach can also 
be applied to calculate the Cronin effect at RHIC [2S]- 

Finally, we point out that even though the mechanism of a hard reaction looks quite differ- 
ently in the dipole formulation from what one is used to in the parton model, the equivalence 
between the dipole approach and the conventional parton model has been demonstrated (for 
single scattering) numerically jTHj and analytically |21| for several processes. 
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